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Abstract 

We review recent work on a continuum, classical theory of thermal fluctuations in 
two dimensional superconductors. A functional integral over a Ginzburg-Landau free 
energy describes the amplitude and phase fluctuations responsible for the crossover 
from Gaussian fluctuations of the superconducting order at high temperatures, to 
the vortex physics of the Kosterlitz-Thouless transition at lower temperatures. Re- 
sults on the structure of this crossover are presented, including new results for 
corrections to the Aslamazov-Larkin fluctuation conductivity. 
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1 Introduction 



The subject of fluctuation conductivity of superconductors is an old one, and 
many classic results can be found in a recent review article by Larkin and Var- 
lamov [1]. Essentially all of this work is restricted to the first corrections to the 
BCS theory, and involving Gaussian fluctuations in the Cooper pair propaga- 
tor. Such an approach is valid provided the Ginzburg fluctuation parameter 
is small. 

In two dimensions, we know that as the superconducting critical temperature 
(Tc) is approached from above, there is ultimately a crossover from such a 
Gaussian fluctuation regime to the vortex physics of the Kosterlitz-Thouless 
(KT) transition. While universal characteristics of the conductivity in the 
latter regime are well understood [2], there is little quantitative understanding 
of precisely how such a crossover occurs. In the context of the small coherence 
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length cuprate supercondTictors, where there is a wide regime of temperatures 
where the Ginzburg fluctuation parameter is large, a quantitative analysis of 
such a regime is surely important in deciphering the importance of pairing 
fluctuations. 

This paper will address the Gaussian-to-vortex crossover in the fluctuation 
conductivity using an approach suggested in Refs. [3,4]; the same approach 
will also be extended to describe the spin-wave-to- vortex crossover at low tem- 
peratures. Using the proximity of the cuprates to a superconductor-insulator 
quantum transition, it was argued that the fluctuation regime could be de- 
scribed by a universal continuum limit of the Ginzburg-Landau free energy. 
This is described by the following partition function for the complex pairing 
order parameter ^(r): 



Zgl = I P*(r)e-^°^/(*^«^) (1) 



2m* ' " ' " 2' 

Here m*, a{T) and h are parameters dependent upon the physics of the un- 
derlying electrons which have been integrated out; the function a(T) vanishes 
at the mean-field transition temperature a{T^^) = 0. Note that the previous 
work [1] was at the level of a mean field treatment of J-'gl, or its one loop fiuctu- 
ations. Furthermore, the strong fiuctuations near the KT transition have been 
invariably treated in previous work in terms of theory for the vortices alone 
or in terms of a phase-only XY model. Refs. [3,4] proposed that the combined 
amplitude and phase fiuctuations implied by (1) should be taken seriously, 
from temperatures well above Tc, and across the vortex physics of the KT 
transition. This paper will present a few results on the electrical conductivity 
associated with the continuum theory of equal-time fiuctuations described by 
(1). Our results for the structure of the crossover contains non-perturbative 
information in the form of a few universal numerical constants which were 
computed in recent computer simulations [3,5,6]. 



2 Review of static theory 



This section reviews a few important static properties of Zql [3,7]. 

A key first point is that the theory (1) cannot predict actual T^. on its own, as 
its value is dependent upon a short distance cutoff. We treat as an input 
parameter, and then find that all subsequent cutoff dependence can be safely 
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neglected. Knowing Tc, it is useful to define the dimensionless coupling g 

\a{T) a{T^y 



9 = 



m*h 



(2) 



which should be a monotonically increasing function of T crossing zero at 
exactly T — Tc. If we take a{T) to be a simple hnearly increasing function of 
temperature {a{T) = ao{T - T^^)), then 



9 = 



1-^ 
T 



(3) 



where Ps{T) is the helicity modulus^ of Zgl (the London penetration depth 
is A^^ = Ane*^ ps{T) / (h'^c^) , where e* = 2e, the Cooper pair charge). In this 
framework, the values of Ps(0) and Tc are the only two needed inputs to predict 
most physical properties of Zql (including the full T dependence of Ps{T)). 
Note that Ps(0) and Tc are, in general, independent of each other, and the 
Nelson-Kosterlitz relation [8] only constraints Ps{Tc)/Tc = 2/n. 

It was shown in Ref. [3,7] that the structure of expansions away from simple 
limiting mean-field regimes of Zgl are best understood in terms of two new 
dimensionless couplings, Q and Qd. These couplings are most useful in high 
and low T limits, respectively. 



For T > Tc, we work with the coupling ^ > given by 
6 1 , / CT 



g^- + -\ni-—), (4) 



where C is a universal numerical constant. Its value is determined by the 

condition that at the critical point with g = 0, Q = Qc = 96.9 ± 3 a universal 
critical value that has been determined in numerical simulations [3,5,6]. Using 
these results, we obtain^ 

C = 79.8 ±2.4. (5) 



A key fact behind our main results here is that for large g (when Q is small), 
all physical properties of Zql can be expanded in a series that involves only 

^ Of course, we do not expect the classical theory in (1) to be valid as T ^ 
because quantum effects will become important; by Ps(0) = —h'^a{0)/{m*b) we 
mean the energy obtained by extrapolating the present approximation for a(T) to 
T = 0. 

^ C is related to the number computed in Refs. [5,6] by C = 6^^^. 
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integer powers of Q. Re-expressing any such series in terms of l/^i by (4) then 
imphes a number of logarithmic terms in g: the coefficients and arguments of 
such logarithms are therefore entirely specified by (4). We should also note 
here that, as discussed in Ref. [7], the combination of (4) and (3) imphes that 
the continuum theory Zql breaks down at a T so large {T > 7rps(0) //cb for the 
T dependence assumed in (3)) that Q ceases to be a monotonically decreasing 
function of increasing T. 

The equation (4) also defines ^ for T < Tc, but its value becomes exponentially 
large as T — > 0. Instead, for T < the useful couphng turns out to be > 
defined by 

g- ^ + ^ ln( ^''^ ] (6) 
Qd tt \GdTcJ 



The coupling Qd vanishes linearly as T ^ 0, and is an increasing function of 
T. At low T we can expand physical properties of Zql in ^d, and again, this 
series involves only integer powers of Qd. The relationship (6) then implies 
logarithms in the expansion in l/\g\- 



3 Conductivity 



A description of electrical transport requires an extension of the partition 
function Zql to a dynamic theory. We do this in the simplest 'model A' 
theory, in which the time evolution of ^(r, t) is given by 



where r is a damping parameter which has no singularities at T = Tg, and C 
is a Gaussian random noise source with the correlator 

(C(r', i')C(r, i)) = 2t^bT5 (r - r') b (t - t') (8) 



Wickham and Dorsey [9] have already computed the conductivity in the theory 
(1, 7, 8) as an expansion in powers of b, and many needed results can be 

obtained simply by transforming their results to spatial dimension d = 2 
(they performed an analysis of perturbation theory in powers of (4 — d) for 
different purposes). As the steps in the calculation are essentially identical to 
theirs, no details of the will be presented here. 
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Our results are most easily expressed in terms of a dimensionless measure of 
the damping parameter, which we define as 



m*h 



(9) 



In BCS theory, the damping arises from decay of Cooper pairs into electrons; 
assuming a linear dependence of a(T) on T as in (3), BCS theory predicts [1] 



Of course, as T — 0, the gapless quasiparticle excitations disappear, and then 
a simple damping term as in (7) is not expected to be valid. Our low T results 
for ZcL should therefore be treated with caution. 

The main new claim of this paper is that the frequency {cu) dependent con- 
ductivity can be written as 

a(a;) = -^T(^,^l (11) 



where T is a universal function (the factor of 12 in the second argument is 
for future notational convenience). The result (11) applies for all values of g 
and Q, but for g <^1, a. different scaling form involving Qd is more convenient 
(see Section 3.3). We discuss the structure of this function for the different 
limiting ranges of g in the subsections below. 



3.1 g>l 



At T > Tc, where g ^ 1, there is an expansion for T in integer powers of 
its first argument Q. The results of Wickham and Dorsey [9], when mapped 
to d — 2, immediately yield an expansion of T to order Q^. For the d.c. 
conductivity, these results are 



The first term in (12) is precisely the old result of Aslamazov and Larkin 
[1]: this can be checked by noting that for g 1^ Q ^ Q/ g from (4), and 
by using (3) and the value for Tbcs in (10)) which yields ultimately a = 
e*^T/(64/i(T — Tc)). The numerical constant A involves a number of different 
contributions, and can be written as 
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-p(«6(0) + /«(0)+/f (0)). (13) 

The first term in (13) is a 'mass' renormalization' correction to the leading 
order term in (12), and arises from (2.1) of Ref. [3] with J = 0.014842966 . . .; 
the remaining terms involving the Fs are associated with 3-loop diagrams 
evaluated by Wickham and Dorsey, and are specified in Eqs. (4.26), (4.28), 
(5.2), (5.3), and (5.4) respectively of Ref. [9]. The momentum integrals in all 
the /'s have to be evaluated in c? = 2. This was done numerically, and we 
obtained 



/^(O) =4.17384 X 10"^ 

7^(0) = -6.92307 X 10-^ 

7^(0)= -1.06867 X 10-^ 
7^(0) =-6.41276 X 10-^ 
7^(0) = -2.93446 X 10"^ 

>1 = 2.28769 X 10"^ (14) 



We can combine (4), (11) and (12) to obtain the following corrections to the 
Aslamazov-Larkin result valid for ^ 1: 



a{0) 



l + —ln{13.3gT/T,) 
ng 



Snhg 

+ ^ (ln2(13.3^T/T,) - ln(13.3^T/T, 



+ 



17287r^ 



9' 



+ 0{l/g') 



(15) 



The overall pre-factor is the Aslamazov-Larkin result. The series is (15) one 
of the main new results of this paper. Using the value of g in (3) as input, it 
yields the T dependence of the fluctuation conductivity. Of course, it is more 
accurate to numerically solve (4) for Q and to then use (12): such a procedure 
would account for all logarithms. 

The results for the frequency dependence of the conductivity are much more 
cumbersome, and we quote only the terms to order Q'^: 



487r y2 



(16) 
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This result is obtained by taking the d — > 2 hmit of Eq. (3.7) in Ref. [9]. 
3.2 g^O 



The Kosterhtz-Thouless transition occurs at Q = Qc, and the system continues 
to be described by (11) across the transition. The singularities of the transition 
appear as singularities in the function T as a function of the Q. The functional 
form of these singularities will similar to those in Ref. [2] as a function of T. 
Of course, there are no arbitrary scale factors in the function T near this 
singularity, but these can only be determination by suitable simulations. 



3.3 5<-l 



For T < Tc, it is useful to expand about a saddle point with 7^ 0, and to 
organize the perturbation theory in powers oi Qd defined in (6). We therefore 
rewrite (11) as 



Again, both Ps{T) and T £, can be expanded in a series involving only integer 
powers oiQo- For PsiT), terms in this series were already computed in Ref. [3]: 

kbT Gd 36 



We computed To to one loop order by the methods of Ref. [9] , and found 



4 Conclusions 



This paper has obtained systematic corrections to the Aslamazov-Larkin fluc- 
tuation conductivity of thermally fluctuating superconductors in two dimen- 
sions. The main results are contained in (11) and (12), and take the form of 
an expansion in integer powers of a renormalized dimensionless coupling Q. 
When Q is re-expressed in terms of bare microscopic parameters (such as the 
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temperature), numerous logarithms appear. However, this microscopic rela- 
tionship is known exactly, and consequently, so is the structure of the all the 
logarithms; this relationship involves non-perturbative information in the form 
of universal constants which were computed in recent numerical simulations 
[3,5,6]. The renormalized coupling Q is related to a bare coupling g by (4), 
and an expansion for the conductivity in inverse powers of g is in (15). The 
coupling g is related to parameters in the Ginzburg-Landau free energy by 
(2). Alternatively, assuming a linear dependence of a(T) on T, we can use the 
simpler expression (3), where Ps(0) = —h^a{Q)/{m*h), is the hehcity modulus 
of the Ginzburg-Landau free energy extrapolated to T = 0. 



It is useful to mention here the expected values of the dimensionless parameter 
Ps{0)/{kBTc) appearing in (3). In BCS theory, we can use standard expressions 
for the parameters in (1), obtained for electrons in a simple parabolic band 
[1]: keeping only the linear T dependence near T^, and extrapolating to T = 
we obtain 



BCS 



nkBTc 



(20) 



where Ep is the Fermi energy. The right- hand- side involves the very large 
ratio Ep/Tc, and this ensures a large regime where g ^ 1 above T^. We 
can also compare to the numerical studies of the cuprate superconductors in 
Ref [10]. They defined a dimensionless parameter t], in terms of which we 
can write PsiS^) / {ksT^ = (2/?7)(T*^^/Tc); for their parameter values we have 
Ps{0)/{kBTc) = 6.8, which also ensures a reasonable large regime over which 
the (7^1 expansion is useful [7]. Finally, we note that near a superfluid- 
insulator quantum phase transition, Ps{0)/{kBTc) is the ratio of two energy 
scales which vanish at the quantum critical point, and so is a universal number 
[3]: this universal number has been computed in a simple model for such a 
transition. 

It would be interesting to extend the present results to other transport co- 
efficients, such as the Nernst effect [10,11]. The present quantitative approach 
is likely to be most useful for observables [7] that do not diverge at Tc. 
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